X-ray imaging with x-ray markers that provide adjunct information but preserve image quality

ABSTRACT

A method and an apparatus for estimating a geometric thickness of a breast in mammography/tomosynthesis or in other x-ray procedures, by imaging markers that are in the path of x-rays passing through the imaged object. The markings can be selected to be visible or to be invisible when the composite markings/breast image is viewed in clinical settings. If desired, the contribution of the markers to the image can be removed through further processing. The resulting information can be used determining the geometric thickness of the body being x-rayed and thus setting imaging parameters that are thickness-related, and for other purposes. The method and apparatus also have application in other types of x-ray imaging.

CROSS REFERENCE OF RELATED APPLICATIONS

This application is a continuation of U.S. patent application Ser. No.12/150,539, filed Apr. 28, 2008, and now U.S. Pat. No. 8,768,026; whichis a continuation-in-part of U.S. patent application Ser. No.11/827,909, filed Jul. 13, 2007, and now U.S. Pat. No. 7,616,801; whichis a continuation-in-part of U.S. patent application Ser. No.11/271,050, filed Nov. 11, 2005, and now U.S. Pat. No. 7,577,282; whichis a continuation-in-part of U.S. patent application Ser. No.10/723,486, filed Nov. 26, 2003, and now U.S. Pat. No. 7,831,296; thedisclosures of which are hereby incorporated by reference in theirentireties.

TECHNICAL FIELD

This patent specification relates to x-ray imaging and more specificallyto accurately finding the geometric thickness of an object being x-rayedthrough including markers in the x-ray path and computer-processing thecomposite markers/object x-ray image to identify the markers and to usethe resulting information to calculate geometric object thickness andimprove imaging.

BACKGROUND

Breast cancer and other breast lesions continue to be a significantthreat to women's health. X-ray mammography currently is the most widelyused tool for early detection and diagnosis, and is the modalityapproved by the U.S. Food and Drug Administration to screen for breastcancer in women who do not show symptoms of breast disease. A typicalx-ray mammography system compresses and immobilizes a patient's breaston a breast platform positioned between an x-ray source and an x-rayimager, and takes a projection x-ray image (called here a conventionalmammogram or simply mammogram) using a collimated cone or pyramid beamof x-rays at appropriate factors such as mA (current), kVp (voltage) orkeV (energy), and msec (exposure time). In the United States, typicallytwo views are taken of each breast, one from above (cranial-caudal, orCC, with the image plane generally at a 0° angle to the horizontal andone from the side (mediolateral-oblique, or MLO, with the image plane atan angle of typically around 45° to the horizontal). Different typicalviews may be taken in other countries. The x-ray source typically is anx-ray tube operating at or in the neighborhood of 25-30 kVp, using amolybdenum, rhodium, or tungsten rotating anode with a focal spot ofabout 0.3 to 0.4 mm and, in some cases, 0.1 mm or less. An anti-scattergrid between the breast and the imager can be used to reduce the effectsof x-ray scatter. The breast is compressed to reduce patient motion andalso for reasons such as reducing scatter, separating overlappingstructures in the breast, reducing the x-ray thickness of the imagedbreast and making it more uniform, and providing more uniform x-rayexposure. Traditionally, the imager has been a film/screen unit in whichthe x-rays impinging on the screen generate light that exposes the film.In the last several years, mammography systems using electronic digitalflat panel array receptors have made significant inroads. A Selenia®digital mammography system with such a digital flat panel x-ray receptoror imager is offered by Lorad, a division of the assignee hereof,Hologic, Inc. of Bedford, Mass. See brochure “Lorad Selenia®” DocumentB-BI-SEO US/Intl (May 2006) copyright Hologic 2006. Digital mammographyhas significant advantages and in time may fully supplant film/screensystems. Additional information regarding digital mammography systemsand processes offered by the common assignee can be found at<www.hologic.com>. Digital tomosynthesis also has made advances and theassignee hereof has exhibited breast tomosynthesis systems at tradeshows and has carried out clinical testing. It is a three-dimensionalprocess in which several two-dimensional projection views are acquiredat respective different angles but lower x-ray dose than conventionalmammograms, and are reconstructed into tomosynthesis slice views thatcan be along any desired plane in the breast. For tomosynthesis, thebreast is still immobilized but may be compressed to the same or lesserextent than in conventional mammography. See, e.g., InternationalApplication WO 2006/058160 A2 published under the Patent CooperationTreaty on Jun. 1, 2006 and Patent Application Publication No.2001/0038681 A1, PCT application International Publication No. WO03/020114 A2 published Mar. 13, 2003, U.S. Pat. Nos. 7,142,633,6,885,724, 6,647,092, 6,289,235, 5,051,904, 5,359,637, and 4,496,557,and published patent applications US 2004/0109529 A1, US 2004/0066884A1, US 2005/0105679 A1, US 20050129172A1, and Digital Clinical Reports,Tomosynthesis, GE Brochure 98-5493, November 1998. Reference markers canbe used in x-ray imaging for purposes such as checking the rotationangle and unwanted shift of center of rotation of an x-ray source andreceptor (imager), and fiducial phantoms can be used in 3D angiographyto calibrate for irregular scan geometries. See, e.g., U.S. Pat. Nos.5,051,904, 5,359,637, and 6,289,235, N. Navab, et al., Dynamicgeometrical calibration for 3D cerebral angiography, SPIE Vol. 2708, pp.361 370, and said PCT published application WO 03/020114 A2. Atomosynthesis system specifically for imaging patients' breast isdisclosed in commonly owned U.S. Pat. Nos. 7,123,684 and 7,245,694. Thesame system can be selectively used for mammography and tomography, inthe same or different compressions of the patient's breast.

It is desirable to know the geometric thickness of the immobilizedbreast in both film/screen and digital flat panel x-ray mammography aswell as in tomosynthesis in order to make appropriate setting for theimaging procedure, such as settings for the x-ray tube that control thex-ray beam. Knowing the breast thickness can also help in quantitativeassessments regarding x-ray images of the breast, such as in assessingthe nature and clinical significance of x-ray attenuation properties ofthe breast. It can also be used in tomosynthesis reconstructions such asto determine the required reconstructed field of view or desired displayfield of view. It has been proposed to use encoders to measure thegeometric height of a breast compression paddle and use the result toestimate the geometric thickness of the breast but this processtypically has a relatively high error because of factors such as tiltingand geometric distortions of the paddle as it compresses the breast, andbecause the encoders must be calibrated. The information on the paddleheight is oftentimes stored in the DICOM header associated with themedical image, and this information is used for several purposes,including the calculation of body part radiation exposure. It also hasbeen proposed to derive breast thickness information from measuring thepaddle compression force versus time and to use the results to guidecontrol factors such as kV, mAs, and filter selection, as in commonlyassigned U.S. Pat. No. 7,123,684. The geometric thickness of any otherbody part being x-rayed also can be of interest, and one process forbasing an estimate on non-contact ranging using ultrasound is discussedin U.S. Pat. No. 4,597,094. It has also been proposed to use acalibration phantom with embedded pellets of high-density materialserving as markers and to use the imaged markers in calibrating thesystem, such as once a week. See U.S. Pat. No. 7,142,633 cited above.

Proposals have been made for automated methods of estimating breastdensity from mammograms, including volumetric and areal estimates. Oneexample, using digitized film mammograms, is given in Br J Radiol. 2006May; 79 (941):378-82 16632617. These methods are understood to derivethe volume of fibroglandular tissue in the breast tissue for across-sectional area in a mammogram. The accuracy of these methods isbelieved to be related to the accuracy of measurement of the actualcompressed breast. In addition, information has been published proposingthat breast density is related to the breast cancer risk, See CancerEpidemiology Biomarkers & Prevention, Vol 13, 715-722, May 2004,American Association for Cancer Research.

The patents and other publications identified above, including thebrochure “Lorad Selenia™” and said published application WO 2006/058160(corresponding U.S. patent application Ser. No. 11/791,601), are herebyincorporated by reference in this patent specification.

SUMMARY OF DISCLOSURE

In one non-limiting example, a breast imaging method comprises x-rayimaging a pattern of markers on a compression paddle and a patient'sbreast immobilized between the paddle and an x-ray imaging receptor toform a composite x-ray image in which the x-ray image of markers iscomposited with the x-ray image of the breast, computer-processing thecomposite x-ray image to derive geometric information regarding theimaged markers, computer-processing the geometric information to derivethickness information related to a thickness of the immobilized breastor other information related to the position of the pattern of markersin space or relative to system components, and using the thickness orspatial position information for improving breast imaging. The imagedmarkers may be invisible within the breast outline when viewing thecomposite image in typical clinical settings or, in alternativeembodiments, they may be visible in the composite image. The method canfurther include removing some or essentially all the contribution of themarkers from the composite x-ray image.

As another non-limiting example, an apparatus comprises a compressionpaddle with a pattern of markers that immobilizes a patient's breastagainst a breast platform, an x-ray source on one side of the paddle andplatform and an x-ray imaging receptor on the other side, a processorreceiving image information from the receptor related to a compositex-ray image of the markers and breast in which the markers are notvisible in a clinical setting, said processor computer-processing thecomposite image to derive geometric information regarding the markerstherein and computer-processing the geometric information to derivethickness information related to thickness of the immobilized breast,and said processor using the thickness information to improve imagingthe breast. The processor can further process the composite image tocompletely or essentially remove some or essentially all contribution ofthe imaged markers.

The method and the apparatus can be applied to both mammography andtomosynthesis systems, and the term “mammography/tomosynthesis” is usedin this patent specification to mean any one of a mammography system, atomosynthesis system, and a fusion system that can selectively carry oneor both of mammography and tomosynthesis, including while a patient'sbreast remains immobilized.

In addition to or instead of being on the compression paddle, similarmarkers can be placed on other objects. For example, markers can beplaced on the magnification table that is placed between the normalbreast platform and the breast to increase the distance between thebreast and the imaging plane and thus magnify the breast image. Markerscan be placed on the normal breast platform if it desirable to calculateor confirm the spacing between the breast platform and the image plane.In general, markers can be placed on any object that is between thex-ray tube and the image plane of the x-ray receptor to help calculatethe distance between the object and the image plane or other objects.

A pattern of parallel lines, like the mm marks of a ruler, can be usedas the pattern of markings. However, many other patterns of markings canbe used and may be preferable depending on the goals of the measurement.For example, a pattern can be dots or a distributed collection of shortincrements of lines that will be imaged as very small areas that arewell dispersed in the x-ray image. The distribution can even bepseudo-random, to further help in making the imaged pattern invisible inclinical settings, so long as the pattern has properties that can helpfind it when imaged in the x-ray image; for example the pattern isknown, or has enough sharp edges that can differentiate the patternelements from background breast images through suitable edge detectionof other computer processing. If the markers in a pattern are elongated,they can be at any desired orientation.

While this patent specification discusses, as one alternative, patternsof markings that when imaged are invisible in the composite x-ray imagewhen viewed in typical clinical settings, it also explains analliterative in which the imaged markers are visible in the compositeimage when viewed in clinical settings. In the first alternative, thex-ray attenuation properties of the patterns of markers can have valuesthat are insufficient to make the imaged markers visible in clinicalsettings but sufficient to allow them to be identified through computerprocessing. In the second alternative, the pattern of markers may havesufficient x-ray attenuation to make the imaged markers visible in thecomposite image in typical clinical settings; however image processingcan process the composite image to remove some or essentially allcontribution of the imaged markers. Even when the imaged markers areinvisible in the composite image in typical clinical settings, theircontributions to the composite image can be removed through furthercomputer processing to reduce possible interference with imageprocessing such as CAD (computer aided detection). Another example ofthe patterns of markers arises in systems capable of tomosynthesisimaging. Tomosynthesis examinations acquire of raw, or projectionimages, and subsequently reconstruction them into tomosynthesis sliceimages that can be desired orientations and represent breast slices ofdesired thicknesses. Because the radiological assessment is typicallymade from the reconstructed slice images, patterns of markings andanalysis as described in this patent specification are useful even ifthey are visible in the tomosynthesis projection images. This is becausetheir contrast and visibility will be greatly reduced in thereconstructed slices that encompass the breast volume, because thepatterns are physically located outside the breast. Because theseexamples relax the requirement for complete invisibility, additionaluses become practical. In a typical digital mammography or tomosynthesisexamination, the setting of x-ray exposure techniques are determinedthrough the use of Automatic Exposure Control (AEC) methods, whichtypically use a low dose x-ray pre-pulse. The low doses employed in thepre-pulse means that the visibility of patterns of markings is reducedand so their contrast will need to be increased accordingly. Patterns ofmarkings and methods according to this patent specification can be usedwhen the appearance of the imaged markers is reduced through imageprocessing, or when the imaged markers do not appear in the breastimages of interest because those images are tomosynthesis reconstructedimages representing slices of the breast while the Imaged markings wouldmainly influence reconstructed images (if any) of slices outside thebreast volume.

As yet another non-limiting example, computer software stored in atangible storage medium can control an x-ray system to carry outimplementation of the process described in this patent specification.Information such as breast thickness, paddle height and paddledeformation so obtained can be stored in the DICOM header associatedwith the medical image for subsequent uses.

The process and apparatus are not limited to applications for breastimaging but are applicable to other x-ray procedures as well.

BRIEF DESCRIPTION OF THE DRAWINGS

The features of the subject matter of this patent specification can bemore readily understood from the following detailed description withreference to the accompanying drawings wherein:

FIG. 1 is a schematic sectional view illustrating a geometry useful inestimating a thickness of a patient's breast immobilized between acompression paddle that has a pattern of markers and a breast platform,when imaged with x-rays in accordance with one embodiment of the presentpatent specification;

FIG. 2 illustrates a two-dimensional projection, such as a mammogram ora tomosynthesis projection image, containing a pattern of markers in theform of parallel bars or lines, where the imaged pattern has beenemphasized for illustrative purposes but would be invisible in aclinical mammogram or other projection x-ray image;

FIG. 3 illustrates an effect of shift-and-add correlation processing ofthe pattern of FIG. 2 where the shift distance does not equal the pitch,or separation distance, between the markers in the image;

FIG. 4 illustrates an effect of shift-and-add correlation processing ofthe pattern of FIG. 2 where the shift distance equals the pitch of themarkers in the image;

FIG. 5a illustrates a ruler pattern that can be used to verify anexample of the process disclosed in this patent specification; FIG. 5billustrates a top view of a compression paddle in current public use inmammography systems offered by the assignee hereof, with automaticexposure control (AEC) markers that are not visible in a clinical image;and FIG. 5c illustrates the compression paddle with an added example ofa pattern (a ruler pattern as in FIG. 5a ) for use with an embodiment ofthe disclosed process;

FIGS. 6a and 6b pertain to x-ray images taken with a compression paddleand ruler as in FIG. 5c , compressing a 4.5 cm thick BR50/50 phantom(BR50/50 is material simulating a breast comprised of 50% fatty and 50%glandular tissue, such as sold by Cirs Inc Model 014A MammographyPhototimer Consistency Testing Slabs, available from CIRS in Norfolk,Va., see http://www.cirsinc.com/pdfs/014Acp.pdf) in the case of FIG. 6abut a 4.5 cm thick cadaver breast phantom in the case of FIG. 6b . Ineach of FIGS. 6a and 6b : strip 1 illustrates a raw x-ray image ofpattern and phantom/breast; strip 2 illustrates a pattern-correctedimage obtained by removing the image of the pattern from the image ofstrip 1: strip 3 illustrates a difference image obtained by subtractingthe strip 2 image from the strip 1 image; strip 4 illustrates ashift-and-add image of raw image 1 where the shift is at the correctdistance and thus the imaged pattern of bars appears in the strip 4image as an interference signal that is being detected; and strip 5illustrates a shift-and-add image of image 2 where the shift also is atthe correct distance but no interference signal is visible,demonstrating that the pattern that was present in strip 1 has beenremoved through pattern correction processing;

FIG. 7a is a graph illustrating an example of amplitude variations of aninterference signal as in FIG. 6b (strip 4) vs. estimated compressionpaddle height above a breast plate at normal 1× clinical x-ray dose. Thetrue height of the compression paddle was 4.5 cm, and the measuredheight using the methods disclosed in this patent was 4.49 cm,identified through the peak of the curve in FIG. 7 a. The interferencesignal amplitude peaks at a point within 0.01 cm of the paddle height.FIG. 7b illustrates the interference signal from a specificshift-and-add image in the same example, as an averaged signal plot inthe spatial domain in the upper plot and in the frequency domain in thelower plot, after frequency analysis of the spatial domain plot throughFFT (Fast Fourier Transform) analysis;

FIG. 8a is a graph Illustrating an example of amplitude variations of aninterference signal as in FIG. 6b (strip 4) but with a 10.0 cm thickcadaver phantom vs. compression paddle height at 1/15^(th) of a normal1× clinical x-ray dose. The true height of the compression paddle was10.16 cm, and the measured height using the methods disclosed in thispatent was 10.13 cm. The interference signal amplitude peaks at a pointwithin 0.2 cm of phantom thickness. FIG. 8b illustrates the interferencesignal from a specific shift-and-add image in the same example, as anaveraged signal plot in the spatial domain in the upper plot and in thefrequency domain in the lower plot, after frequency analysis of thespatial domain plot through FFT (Fast Fourier Transform) analysis;

FIG. 9 is a partial side view of a mammography system in which apatient's breast is immobilized between a compression paddle and abreast platform and is imaged on a flat panel digital imager, with orwithout the use of an anti-scatter grid, in accordance with oneembodiment of the system and method disclosed in this patentspecification; and

FIG. 10 is a block diagram illustrating main portions of a system inwhich the disclosed process is used.

FIG. 11 illustrates an example using two sets of markers, at differentheights above an imaging plane (detection layer) that can facilitatecalculation of the heights, and illustrates an x-ray image showing theimaged markers (that would not normally be visible in a clinical image);

FIG. 12 is a flowchart illustrating steps in a process of estimatingheight of objects with markers;

FIG. 13 illustrates shift-and-add according to the flowchart of FIG. 12;

FIG. 14 illustrates two different patterns of markings: a pattern “a” inthe form of a regular pattern of parallel short lines and a pattern “b”that can be arbitrary and can even be random or pseudo-random, so longas it has characteristics that can differentiate it from a backgroundx-ray image to which it contributes; and

FIG. 15 is a schematic Illustration of another way of measuring bothheight and tilt of a compression paddle.

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS

This patent specification describes methods and systems in which thegeometric thickness of a body part such as the breast that is beingx-rayed is accurately determined and used to improve imaging in a mannerthat does not inconvenience the patient or the health professional. Italso describes similar processes to determine other thicknesses,heights, or distances between objects in x-ray examination procedures.

In describing preferred embodiments illustrated in the drawings,specific terminology is employed for the sake of clarity. However, thispatent specification is not intended to be limited to the specificterminology so selected and it is to be understood that each specificelement includes all technical equivalents that operate in a similarmanner. In addition, a detailed description of known functions andconfigurations will be omitted when it may obscure the subject matter ofthe invention described in the appended claims.

As illustrated in the example of FIG. 1, a mammography or tomosynthesissystem immobilizes a patient's breast 12 between a compression paddle 10and a breast platform 14. An x-ray source 16 when energized emits x-raysfrom a focal spot 18 that pass through paddle 10, breast 12 and platform14 and are imaged at an imaging plane of an imaging receptor 20 that isbelow platform 14 and may for some x-ray procedures be separatedtherefrom by an anti-scatter grid not shown in FIG. 1. A pattern 22 offive markers “a, b, c, d, e” is formed in or on paddle 10 and is imagedby the x-rays as respective imaged markers “A, 8, C, D, E” in theprojection breast image formed at the imaging plane of receptor 20. Theoverall geometry of a mammography/tomosynthesis system to which thearrangement of FIG. 1 can be added is described in said publishedapplication WO 2006/058160.

It can be appreciated that if the geometry of markers 22 is known andthe geometry of imaged markers 24 is measured in the x-ray image,trigonometric calculations can yield the distance between markers 22 andthe imaging plane. The distance between markers 22 and the bottom ofcompression paddle 10, and the distance between the imaging plane andthe top of breast platform are known. These distances are sufficient tocalculate the geometric thickness of immobilized breast 12. Anadditional challenge arises in cases where it is not acceptable ordesirable to have an easily visible contribution of the imaged markersto the composite x-ray image. For example, it may not be practical in atleast some cases to have imaged markers visible in the breast image asthey may detract from the clinical value of the x-ray image.

The new approach disclosed in this patent specification involves, as onealliterative, forming a composite x-ray image of tissue and superimposedmarkers such that the markers are invisible in the composite image in aclinical setting but can be identified though computer-processing of thecomposite image. As another alternative, the markers are visible in thecomposite image. In either or both alternatives, the disclosed processand system further provide a way to remove contributions of the imagedmarkers from the composite image, if desired. “Invisible” in this patentspecification means that the presence of imaged markers are not visibleto an observer within the outline of a patient's breast in a realisticanatomical breast background under clinical x-ray dose, in typicalclinical settings and in the type of visual examination of breast x-rayimages that is done in conventional clinical screening and/or diagnosis.The approach disclosed in this patent specification is applicable toconventional mammography using flat panel or other digital imagingreceptors, to digitized film, to tomosynthesis imaging using suchreceptors, and to other types of x-ray procedures and other x-rayimaging receptors, where it is desirable or helpful to accuratelyestimate the geometric thickness of the body or body part being x-rayedand to use the results to improve imaging or for other purposes.

A non-limiting example of the new approach uses one or more patterns ofmarkers 22 that produce imaged markers 24 whose contrast relative to thebreast image, or image of other body parties, can be so low that theimaged markers would not be seen in conventional assessment of x-rayimages but can be detected through computer-processing of the imagesusing techniques such as correlation analysis. This is made possible byusing markers 22 that are in patterns with appropriate properties suchas spatial distribution or other spatial and x-ray density propertiesas, for example, repeated line bars at known pitch. Correlation analysisof the imaged markers 24 then can allow detection of the position andpitch of imaged markers 24 in the composite markers/breast image. Oncethe geometry of the imaged markers 24 has been found, further computerprocessing can remove the contributions of markers 24 from the compositeimage. This removal can be so effective that an image that should beflat is in fact flat despite the presence of markers 22 at paddle 10.“Flat image” here means an x-ray Image taken with a uniform phantombetween paddle 10 and image receptor 20, i.e., and image that shouldhave a uniform gray level within system tolerances and when viewed intypical clinical setting appears to be flat.

The effectiveness and robustness of the new approach can be demonstratedby experiments with a common ruler made of transparent plastic andhaving small painted black bars marking the millimeters (mm). In theseexperiments, the ruler is placed on paddle 10, and the ruler's mm marksserve as a pattern of markers 22. Breast phantoms of geometric thicknessfrom 0 cm to 10 cm are immobilized between the paddle and the breastplatform and are imaged under clinical x-ray dose levels in amammography unit such as the system available in this country under thetradename Selenia™ from the assignee hereof, Hologic, Inc. of Bedford,Mass. Such experiments can demonstrate that geometric breast thicknesscan be calculated with errors within a ±1 mm range. Under clinical x-raydose and realistic anatomical breast background, the imaged mm marks ofthe ruler are invisible to an observer. In a flat field image with auniform phantom, the ruler image is faint but perceptible. However,additional computer processing as disclosed in this patent specificationcan remove this residue of the ruler image to thereby leave no visibleartifacts in both flat field images and breast images. In similar testswith 1/10th to 1/20th of normal clinical x-ray dose, the ruler marks canstill be identified through computer processing and breast thickness canbe measured to within ±2 mm error range. This can allow geometric breastthickness to be computed from an x-ray image obtained with the AEC(automatic exposure control) scout shot, the very low dose x-ray shottypically taken for automatic exposure control purposes in mammographyunits such as the aforementioned Selenia™ system before actual beastimaging. The breast thickness calculated from the AEC shot informationcan then be used for more accurate AEC procedures to achieve betterquality breast images than with less accurately determined geometricbreast thickness, and/or for other purposes.

In the illustrative examples of FIGS. 1-4, five tiny bar markers 22 areused at paddle 10, uniformly spaced at a known pitch “d.” They projectin the image formed at receptor 20 as five similarly shaped imagedmarkers 24 that also are uniformly spaced but at an unknown pitch “D”.As illustrated in FIG. 1, “h” is the distance between the focal spot 18in x-ray source 16 and compression paddle 10, “H” is the distance fromfocal spot 18 and the image plane of receptor 20, the breast thicknessis designated “breast_thickness,” and “s” designates the distancebetween the breast platform and the image plane. Parameters “d” and “H”and “s” are known, and parameter “D” can be determined as explained inthis patent specification. Parameter “h” then can be calculated from therelationship h=(d/D)H, and breast thickness can be calculated as(breast_thickness=H−h−s). In this example, the centerline of the x-raybeam is normal to an centered at the image plane 20; however, if thecenterline is at a different angle and point of incidence in the imagingplane, as in the case of some tomosynthesis projection images, a uniformruler pattern as in FIG. 1 may be geometrically distorted. However, thedistortion can be calculated at least approximately from knowing systemgeometry at the time the projection image is taken, e.g. from variousencoders for the position of the x-ray source and from the constructionof the system, so these geometric distortion can be effectivelyaccounted for in the calculations disclosed in the patent specification.

In this alternative example imaged markers 24 would not be visible in arealistic breast image in a clinical setting, but to assist in thediscussion below they are illustrated as visible in a compositemarker/breast image in FIG. 2. Correlation processing can be carried outby shifting the image of FIG. 2 by a specified distance, e.g. to theright, adding the original and shifted images, then shifting the imageof FIG. 2 again by the same distance and adding each shifted image tothe result of the previous shift-and-add, etc., as long and the originaland the last shifted images partially overlap. The shifting in generalis not by an integer number of pixels so appropriate image interpolationcan be carried out to generate the final shift-and-added image. Theresult of this shift-and-add process, or a normalized version thereof,can be expressed as amplitude related, for example, to an averageattenuation or gray level of pixels in the summed images. Thisshift-and-add process is repeated for each of a number of differentshift distances, e.g., for each of a progression of distances thatdiffer in small steps and are in a reasonable range related to expectedbreast thickness ranges and taking into account possible magnificationimaging in which the breast is spaced from platform 14 by a selecteddistance and, if desired, possible geometric distortions of pattern 22because of factors such as a non-normal angle between the x-rays thatform the image of a marker and the imaging plane. If the shift distancematches the actual pitch “D” of the imaged markers 24, the amplitudewould be at a peak because the imaged markers 24 in effect line up andadd while the rest of the composite x-ray image is more random and doesnot increase in amplitude as much. If the shift distance does not matchthe pitch “D,” the amplitude resulting from addition is less becausemarkers 24 are not aligned or at least not as aligned.

The shift-and-add process has been described qualitatively above, but inactual practice it can be performed by processing the pixel values ofthe composite markers/breast image. Consider the pixel values derivedfrom receptor 20 representing pixels that are all within the middleimaged marker 24. After a single shift-and-add by “D,” the pixel valuescontributed by the same marker 24 will be doubled; after anothershift-and-add they will be tripled; and so on. However, if the shift isby a distance different from “D,” at least some of the pixel values forthe marker will not be doubled—they will be added to pixel valuesrepresenting only the breast image. The pixel values for the entirecomposite image can be added or averaged, and the resulting amplitudeshould be highest for shift-and-add processes where the shift distancewas “D.” Preferably, the shift-and-add process is applied to one or moreregions of interest (ROI) likely to contain imaged markers 24 ratherthan to the entire x-ray image. Many mathematical techniques are knownfor use in searching for weak but periodic signals, e.g. signals thatare below a noise level of their environment. Such techniques generallyinvolve some form of correlation analysis, and many correlationprocesses known in mathematics may be used in place of the shift-and-addexample discussed in more detail in this patent specification. However,to the inventors' knowledge no such techniques have been applied in aprocess and apparatus as disclosed in this patent specification, e.g.,for accurate determination of the geometric thickness of breast or othertissue being x-rayed and to use of the resulting information accordingto this patent specification.

A single pattern of markers can be used, or several patterns can be usedat the same level or at different levels. While in these experiments aruler mm marks are used as an example of a pattern, other patterns canbe used so long an they have known geometries and are correlated, e.g.,repeating at known pitch or otherwise correlated. While FIG. 1illustrates a pattern 22 of only five bars, typically many more bars orother marks would be used. Only a selected part of the composite imagecan be processed that is likely to contain a significant portion of thepattern or patterns. For example, if one or more patterns are placed atthe central portion of paddle 10, say 1 cm from the paddle edge that ispressed against the patient's chest wall, then only a strip of the x-rayimage that is likely to contain an image of that portion of paddle 10can be computer-processed as described above. In addition, in caseswhere paddle tilting is deliberately provided, patterns or sub-patternsof markers 22 can be provided at different places of paddle 10, and thethickness of the immobilized breast can be calculated for the portion ofthe breast under each pattern and interpolated for other portions of thebreast. The resulting information can be used to set exposure parametersand possibly for breast image processing purposes specific to respectiveportions of the imaged breast.

FIG. 3 shows four shifts of the image of FIG. 2 where each shift is by aconstant distance that is not equal to the pitch “D,” and FIG. 4 is asimilar illustration but for shifts by a distance “D.” As seen, theimaged markers 24 are misaligned in FIG. 3 but not in FIG. 4, where thesecond marker from the left is the sum of marker 2 of the original imageand marker 1 of the image after one shift by “D,” etc.

An illustrative but non-limiting way to implement the process inpractice involves the following basic steps:

-   -   a. Compression paddle 10 is painted with markers 22 in a line        bar pattern with a known bar pitch, e.g., 1 mm, as in a normal        plastic ruler. Each line bar is about 6 mm long and is extends        in a direction normal to the chest wall such that the bars are        repeated along the chest wall direction in a mammographic or        tomosynthesis procedure. Each paint line is thin, e.g., about        0.1 mm to 0.2 mm, for better correlation analysis. The contrast        of the paint lines preferably should be low enough such that the        imaged markers 24 in the composite marker/breast image are        invisible when the composite x-ray image is viewed under typical        clinical conditions;    -   b. Because the height of paddle 10 can be known to within about        ±1 cm from conventional encoder information supplied in digital        form to the system computer in a mammography unit such as the        Selenia™ system, it can be calculated through basic        trigonometric relationships approximately where the imaged        markers 24 will be in the composite marker/breast x-ray image        and approximately what the pitch will be of imaged markers 24.        These calculations identify a region of interest (ROI) in the        composite image that will likely contain imaged markers 24 or at        least a good portion of the markers, and also identify a likely        shift distance for the correlation analysis or at least a range        of likely shift distances. In addition, if angular encoding is        used to measure tilt of the paddle (see FIG. 15 and description        thereof below) and several sub-patterns are used (as in FIG. 5c        ), the tilt information can be used to further refine the        estimate of the likely location in the image of imaged markers        of each sub-pattern, using known trigonometric relationships;    -   c. In one example of correlation image analysis, this ROI is        shifted incrementally along the chest wall direction that is        normal to line bars, with different shift amounts within the        expected range, and the image resulting from the additions is        checked for the line bar interference pattern;    -   d. Once the shift amount is equal to the line bar pitch, the        resulting image will show strongly enhanced line bar patterns,        say a factor of 50-100 times stronger relative to the        background. In addition to shift-and-add along chest wall        direction, all pixels along the 6 mm long line bar as imaged in        the composite image are added. This can give another factor of        100-300 times more signal in the final line profile;    -   e. In frequency domain analysis, the strongest interference        signal can be picked up at the expected frequency, and the        corresponding shift amount can be found that is equal to the        line bar spacing. From the measured pitch “D” of the imaged        ruler markers at unknown height and the known pitch “d” of the        ruler, the system magnification ratio can be calculated, and        then the paddle height or breast thickness is derived at the ROI        region. It multiple clusters or sub-patterns of line bars are        placed at paddle 10, the analysis can be repeated to get the        paddle height profile and breast thickness at each ROI within        the whole detector area, and to calculate paddle tilt angle or        regional deformations as well;    -   f. The process incorporates principles similar to those used in        the design of a “phase lock amplifier” in analog electronics to        identify a weak electronic signal of known frequency in the time        domain. Similar principles can be used in the spatial domain as        a form of a “phase lock amplifier” to measure a spatial        frequency in an x-ray image;    -   g. Once the pitch of imaged markers 24 and the paddle height are        calculated, the location and the signal magnitude of imaged        markers 24 in the composite marker/breast image can also be        found through conventional trigonometric relationships.        Additional artifact correction algorithms can therefore be        applied to remove the residue signal of the ruler marks (imaged        markers 24) from the composite image, which will make the ruler        marks become completely or essentially invisible even in the        flat field image of the ruler. After such removal of imaged        marker contributions, the x-ray images can be more suitable for        processing such as CAD (computer aided detection) analysis.

FIG. 5a shows a magnified ruler pattern that can be used in theexperiments referred to above. As seen, the ruler pattern has mm marks.The pattern was printed, using a laser printer, onto a thin sheet ofplastic, but can be provided by other printing techniques, or can beengraved into the paddle, or can be formed in some other way, so long asit has the required x-ray properties that would cause imaging thepattern into a composite pattern/breast image as an imaged pattern thatis invisible within the breast outline when viewing the composite imagein clinical settings or, in alternative embodiments, the imaged patternis visible if the composite image is viewed in clinical settings.Optionally in either alliterative, contributions of imaged markers canbe removed from the composite x-ray image through further computerprocessing. FIG. 5b shows a typical compression paddle 10 that can beused with a Selenia™ mammography system, with painted marks related toAEC features, which AEC marks are not visible in the breast x-ray imagein clinical settings though they can be clearly seen on the paddle. FIG.5c illustrates the paddle of FIG. 5b with an attached plastic sheetcontaining ruler patterns similar in principle to that of FIG. 5a butthe pattern in FIG. 5c is in the form of nine different sub-patterns ofruler marks, where both the AEC marks and the imaged ruler marks of thesub-patterns are visible on the paddle but can be invisible when imagedin a breast x-ray image in a clinical setting. For clarity, the marks ofthe nine sub-patterns shown in FIG. 5c are coarser and the parallellines are fewer in number in each sub-pattern than they typically wouldbe in actual practice. While the sub-patterns are shown as extendingonly in the lateral direction, in practice the sub-patterns may extendin different directions that can differ between sub-patterns. And, thesub-patterns need not be in the form of ruler marks, as discussed inmore detail below and different sub-patterns can have different marks indifferent arrangements.

FIG. 6a illustrates results of experiments using a paddle 10 with asuperimposed ruler pattern as seen in FIG. 5c immobilizing a 4.5 cmBR50/50 breast phantom, and FIG. 6b is similar but illustrates resultswhen using a 4.5 cm cadaver breast phantom, in each case with the use ofa Selenia™ mammography system. Strip 1 of FIG. 6a Illustrates an x-rayimage (actually a strip taken from an image) that contains the combinedcontribution of the breast phantom and the ruler mm marks serving as apattern of markers 22. The imaged markers 24 are not seen in strip 1 ina typical observation in clinical settings. Strip 2 illustrates the sameimage but with imaged markers 24 removed using the process disclosed inthis patent specification. No difference can be seen between strips 1and 2 in typical clinical setting observations. Strip 3 illustrates adifference image obtained by subtracting the images of strips 1 and 2from each other. As expected, the imaged markers 24 can be seen in thisdifference image, demonstrating that there is information about markers24 in the image of strip 1. Strip 4 illustrates the result of ashift-and-add operation on the image that contains strip 1, using ashift distance that matches the actual pitch of imaged markers 24. Theshift-and-add operation reinforces and emphasizes imaged markers 24,making them visible in the composite x-ray image and thus in another waydemonstrating that they are present in the composite image that containsstrip 1. In practice, this process detects the interference signal dueto the presence of markers 24 in the composite image of strip 1. Strip 5illustrates the result of a shift-and-add operation on the image ofstrip 2, carried out in the same way as the one that produced strip 4.However, imaged markers 24 are not seen in strip 5, demonstrating thatimage 2 was in fact pattern-corrected by removing essentially all thecontributions of the imaged marker 24 from the composite image of strip1. FIG. 6b has corresponding respective strips 1-5 but pertains toimaging a cadaver breast phantom.

FIG. 7a shows a plot of an amplitude of the interference pattern ofimaged markers 24 (seen in FIG. 6a , strip 3) versus different shiftdistances and also shows the measured breast thickness values, in thecase of using a 4.5 cm thick cadaver breast phantom and normal x-rayimaging dose for mammography. Each of the data points is for a differentshift distance. The peak signal amplitude corresponds to calculatedbreast thickness. As seen in FIG. 7a , the calculated breast phantomthickness is within 0.01 cm of the 4.5 cm phantom thickness. The upperplot in FIG. 7b illustrates for the geometry of FIG. 7a an averagedinterference signal in the spatial domain of imaged markers 24 from onespecific shift-and-add image shown in strip 4 of FIG. 6b . The lowerplot in FIG. 7b shows results from FFT (Fast Fourier Transform)analysis, including a peak signal at the expected spatial frequency.FIGS. 8a and 8B are similar plots, respectively, but for the case ofusing a 10.0 cm thick cadaver phantom and 1/15^(th) of normal x-ray dosefor mammography. The amplitude of the interference signal(peak-to-background) is reduced about 40 times relative to the case ofFIG. 7a (from about 0.2 to about 0.005) but the process disclosed inthis patent specification can still detect this weak pattern signal andcalculate breast phantom thickness to within 0.2 cm of the 10.0 cmphantom thickness, demonstrating the high robustness of the process.

FIG. 9 illustrates certain components of a system that can be used formammography or tomosynthesis, in each case using the process ofdetermining breast thickness and using the result to improve breastimaging. An x-ray source 1 is at one end of a generally C-shaped frame 7and a flat panel digital x-ray imaging receptor 5 is at the other end.X-ray source 1 includes a collimator schematically illustrated at 40 toconfine an x-ray beam 30 emitted from source 1 to a desired footprint atreceptor 5, typically no larger than the area of receptor 5 andpreferably just enough to image a patients breast 3 or at least aselected part thereof, as compressed toward receptor 5 by a compressionpaddle 2 mounted on an arm 6 that in turn mounts to frame 7. A lowerplatform 11, often called a breast tray or platform, is immediatelybelow the breast, and a scatter-reducing grid 4 is between breast tray11 and x-ray receptor 5 and is housed in the same enclosure 12 with thereceptor. As is known in the art, frame 7 can rotate between horizontaland vertical directions of x-ray beam 30. In use for a CC view, paddle 2and its supporting arm 6 are moved up, breast 3 is positioned on tray 11and compressed by bringing paddle 2 down as needed. With suitablecollimation by collimators 40 (which typically collimate in twodirections, of which only one is illustrated in FIG. 9), beam 30 fromsource 1 images the breast onto receptor 5 and the resulting electronicimage information is transmitted to a viewing station 112 (FIG. 10). Theimage typically is rectangular. Preferably, the collimation is such thatbeam 30 illuminates an area of receptor 5 just large enough to show theimage of breast 3, or at least a selected part thereof. Different sizesand shapes of paddles 2 can be mounted to arm 6, and the paddle can beselectively positioned off-center relative to proximal edge 5 a ofreceptor 5 (the left edge in FIG. 9). For mammography, grid 4 can be inthe position shown in FIG. 9. For tomosynthesis data acquisition, grid 4can be retracted to the right of FIG. 9 and out the path of x-rays 30,and a number of projection breast images can be taken at differentangles of the x-ray beam relative to the breast and processed to derivetomosynthesis slice images. Paddle 2 of FIG. 9 can be provided with oneor more patterns or sub-patterns of markers 22 as with paddle 10 ofFIG. 1. Breast platform 11 and receptor 5 of FIG. 9 correspond toplatform 14 and receptor 20 of FIG. 1.

A system that embodies an example of the new approach to finding andusing geometric breast thickness is illustrated in block diagram form inFIG. 10. The system can be a mammography system, or a tomosynthesissystem, or a fusion system that can be used for either or both purposes.A beast immobilizing device 102 is between an x-ray source 104 and animage receptor/grid unit 106. Device 102 can comprise paddle 10 andbreast platform 14 of FIG. 1, or paddle 2 and breast platform 11 of FIG.9. Unit 104 can comprise source 16 and focal spot 18 of FIG. 1, orsource 1 and collimator 40 of FIG. 9. Unit 106 can comprise receptor 20of FIG. 1 or the combination of receptor 5 and grid 4 or FIG. 9. Acontrol unit 108 is coupled with units 102, 104 and 106 to control theiroperation and to send and receive information regarding their operation.An image processor 110 is coupled with unit 106 to receive image datatherefrom and exchange control and other information therewith, and alsois coupled with unit 108 to exchange control and other information. Adisplay/operator console unit 112 is coupled with image processor 110and control 108 to exchange display and other information therewith, andto display images provided by unit 112 and receive user input throughappropriate interface devices. The units illustrated in FIG. 10 canInclude the functionalities of the corresponding units of thecommercially available Selenia™ mammography system as well as theadditional processing functionality to carry out the process ofdetermining breast thickness and using the result for AEC and/or otherpurposes. The additional processing functionality can further includebreast density calculations that use information regarding the geometricthickness of the breast as imaged at one or more regions of the imagingreceptor.

An initial application of the disclosed method and apparatus is toprovide accurate estimates of breast thickness in mammography and/ortomosynthesis systems using the AEC scout view before the main exposureand to use the thickness information to set x-ray technique, to provideaccurate estimates of breast thickness for more accurate breast densitycalculations, and for other purposes such as control over tomosynthesisslice image reconstruction and display. In tomosynthesis, the thicknessinformation from a projection image can be used in AEC for otherprojection images and/or in AEC for a subsequent mammogram taken in thesame breast compression in the same system. Similarly, the thicknessinformation from a mammogram can be used in AEC for taking subsequenttomosynthesis views or other x-ray images taken in the same breastcompression in the same system. Another use of the process and apparatusfor compositing invisible or visible markers into an x-ray image and ifdesired removing such imaged markers through further processing is indigital encoding of selected information into an x-ray image. Thecomposited marker information can be used for a wide variety of purposessuch as identifying the patient, the procedure for acquiring the datafor the image, encoding other information regarding the image, or forproviding any other type of information embedded by way of invisiblemarkers. A similar process and apparatus can be used in other x-rayprocedures to determine the geometric location or thickness of the bodyor body part being x-rayed, as in bone densitometry and generalradiography. Still another use of the disclosed process and apparatus isfor registering mammograms or other x-ray images, where the imagedmarkers that can be visible or invisible, as appropriate for aparticular application, can be used as registration marks that can beleft in the registered images when they are invisible in normalobservation and this is desirable, or can be left as visible marks, orcan be removed in either case through further processing as describedabove. As one non-limiting example, the breast thickness calculated asdescribed in this patent specification can be displayed separately fromor together with the breast x-ray images, e.g., can be displayed to thetechnician operating a system such as Selenia™ on a numerical display.

Other examples of applications of the new approach include using two ofmore patterns of markers, for example using two patterns or sub-patternsas illustrated in FIG. 11, where one pattern is above the patient'sbreast and another is below the breast, such that the two patterns donot overlap or only partly overlap in the composite x-ray image. Theupper pattern can be on the otherwise conventional compression paddle ofa mammography/tomosynthesis system and the lower pattern can be on thetop surface of an otherwise conventional magnification table that isplaced between the breast and the breast platform when it is desired tomagnify the breast x-ray image. FIG. 11 includes an illustration to theright of a composite x-ray image that would result from using twoexamples of such patterns. The markers in one alliterative embodimentwould not be visible when viewing the image in a clinical setting. Inanother alternative embodiment, the imaged markers can be visible (whenappropriate patterns 22 are used) and can be removed when it is desiredto view the composite x-ray image without visible residue of imagedmarkers.

FIG. 12 is a flowchart Illustrating one example of a process accordingto the new approach disclosed in this patent specification, and showsthe following main steps of a process for calculating the thickness ofan object that is between an upper and a lower plane and is beingx-rayed. At least one pattern of markings is present at least at one ofthe upper and lower planes. In a non-limiting example, the upper planecan be at a compression paddle and the lower surface at a breastplatform in a mammography/tomosynthesis system. The main steps can be:

-   -   a. Acquire the x-ray image, with a pattern of markers imaged        within the field of view, e.g., as illustrated in FIG. 1 or 11;    -   b. Identify the approximate region of interest (ROI) where the        imaged markers are approximately expected in the x-ray image.        This can be derived through conventional calculations from        knowing the locations of the pattern(s) relative to the        acquisition geometry, including from a conventional encoder        providing information regarding the approximate height of a        compression paddle;    -   c. Start the calculation at an estimated thickness t=t₀ that is        less than the true thickness of the object. If the system        reports an approximate thickness of the object using, for        example, an encoder for paddle height, the encoder value can be        used. For example, if the encoder reports the paddle height is 5        cm, a starting height could be 4 cm. If there is no known        estimate for the true paddle height, the starting height could        be zero;    -   d. Generate a shift-and-added x-ray image. The amount of        shifting per image is determined from the thickness t being        estimated, and from system geometry, e.g., according to the        expressions relating thickness to known parameters discussed        above in this patent specification;    -   e. Generate a profile by summing the image in a direction        perpendicular to the lengths of the imaged markers;    -   f. Calculate the one dimensional FFT of the profile;    -   g. Find the amplitude of the peak;    -   h. Increment to the next thickness to test. Typical increments        depend upon the desired final thickness measurement accuracy,        and for digital mammography/tomosynthesis can be, as a        non-limiting example, in the range of about 0.1 mm to 1 mm;    -   i. If the thickness is less than the desired maximum thickness        to test, repeat at step d;    -   j. Now find the thickness that generated the largest peak        amplitude in step g. This thickness is the calculated thickness        of the object according to the pattern above the object. Similar        processes can be used if the parameter to be calculated is the        position in space or with respect to another portion of the        system, of the upper or lower patterns of markings, or the        difference in position between the two patterns of markings, or        if a parameter is sought that is related to the position in        space or in relation to a system component, or each or two or        more patterns of markings that are approximately at the same        level, such as two or more patterns on a compression paddle that        might tilt or deform in the course of compression.

The process described immediately above can be adapted to the moregeneral case where it is desired to calculate the position in space ofone or more patterns of markings that are between an x-ray source and animage plane of an x-ray imaging receptor.

FIG. 14 illustrates an upper pattern of markings “a” that is similar tothe pattern in FIG. 5a and to a sub-pattern of FIG. 5c , and a lowerpattern of markings “b” consisting of the lines of pattern “a” brokendown into four rows each consisting of much shorter lines. Because anyarbitrary pattern whose general shape of certain other properties areknown can be detected in the x-ray image through spatial domaincorrelation processes (without FFT) as long as the imaged pattern hasenough signal, the arbitrary pattern can be found with sufficientaccuracy according to the approach disclosed in this patentspecification. Even a random or pseudo-random pattern of small areas inthe composite x-ray image can be detected if it has properties thatgenerally distinguish it from the image of a breast or other object inthe composite x-ray image. One example of such characteristics can beedges that are too sharp and appear too frequently to belong to an imageof a breast or another object. In the example of pattern “b” in FIG. 14,The x-ray imaged rows of shorter lines of pattern “b” can be shiftedrelative to each other through computer processing in a directionperpendicular to their lengths to conceptually collapse them intopattern “a.” As long as pattern “b” has approximately the same number ofpixels as pattern “a,” It can be found nearly equally well in thecomposite x-ray image. Pattern “b” need not consist of short lines;other shapes can also work. A pattern such as “b” can ensureinvisibility in the composite x-ray image more reliably and under morerelaxed viewing conditions.

FIG. 15 schematically Illustrates measuring not only the height of acompression paddle such as paddle 2 of FIG. 9 but also tilt of thepaddle as it compresses the patient's breast against breast platform 12.As illustrated, support 6 for paddle 2 is coupled with linear andangular encoders 90 in addition to being supported for vertical motionalong column 7. The linear encoder can be an optical encoder or anencoder as used in the Selenia™ system available from the commonassignee. The angular encoder can be any suitable optical encoder or aresolver that can measure the angle of paddle 2 relevant to a referenceas the paddle compresses the patient's breast and tilts from a planeparallel to breast platform. Both the linear and angular encodersprovide digital information to control 108 and image processor 110,which information can be used as discussed above to assist inidentifying regions of interest to which correlation processing can beapplied. A linear encoder of this type has been in public use in thiscountry for more than a year. An angular encoder for breast compressionpaddle that provides digital information to a control such as 108 hasnot been in public use but a mammography system with a built inpotentiometer that is not providing digital Information and is not usedto measure tilt and provide it for use in the system has been in systemsshown in at least one trade show in this country.

Software that controls the process described above can be stored in atangible computer readable storage medium to be used as a computerprogram product and/or can be transmitted via a computer network orother transmission medium.

The above specific embodiments are illustrative, and many variations canbe introduced on these embodiments without departing from the spirit ofthe disclosure or from the scope of the appended claims. For example,elements and/or features of different examples and illustrativeembodiments may be combined with each other and/or substituted for eachother within the scope of this disclosure and appended claims. As one ofmany examples, the processes disclosed and claimed in this patentspecification can be applied to CR based x-ray imaging, where the chargestorage phosphor or x-ray detection material is not rigidly fixed inrelative to the breast platform, or similarly to digitized film images.

The invention claimed is:
 1. A method comprising: concurrently x-rayimaging a breast compression paddle comprising at least one paddlemarker and a patient's breast compressed between the compression paddleand an x-ray imaging receptor; obtaining an x-ray image in which the atleast one paddle marker and the breast are imaged, wherein an image ofat least a part of the at least one paddle marker is superimposed on animage of the breast in the x-ray image, and wherein the at least onepaddle marker is invisible within an outline of the breast in the x-rayimage; identifying at least one region of interest in the x-ray imagethat corresponds to at least a portion of the at least one paddlemarker; processing the at least one region of interest so as to derivegeometric information regarding the at least one paddle marker image inthe x-ray image; and processing the geometric information to deriveinformation regarding the compression paddle.
 2. The method of claim 1,wherein the information regarding the compression paddle comprises atleast one of a paddle height, a paddle deformation, a paddleorientation, and a compressed breast thickness.
 3. The method of claim2, further comprising storing the information regarding at least one ofthe paddle height, the paddle deformation, the paddle orientation, andthe compressed breast thickness in association with stored imageinformation regarding the imaged breast.
 4. The method of claim 1,wherein the at least one paddle marker comprises a plurality of paddlemarkers.
 5. The method of claim 4, wherein the plurality of paddlemarkers are disposed along a width of the compression paddle.
 6. Themethod of claim 5, wherein the geometric information comprises a paddleheight profile at each of the plurality of paddle markers.
 7. The methodof claim 4, wherein at least some of the plurality of paddle markers arevisible in the x-ray image when viewed in clinical settings.
 8. Themethod of claim 4, further comprising selecting the region of interestin the x-ray image that is likely to contain at least a substantialportion of the plurality of imaged paddle markers.
 9. The method ofclaim 1, wherein the x-ray imaging receptor comprises at least onereceptor marker, and wherein the at least one receptor marker is imagedconcurrently with the breast compression paddle and the patient'scompressed breast.
 10. A method comprising: performing a first x-rayimaging acquisition concurrently on a breast compression paddlecomprising at least one paddle marker and a patient's breast compressedbetween the compression paddle and an x-ray imaging receptor, whereinthe first x-ray imaging acquisition is performed with a first x-raydose; obtaining an x-ray image in which the at least one paddle markerand the breast are imaged, wherein an image of at least a part of the atleast one paddle marker is superimposed on an image of the breast in thex-ray image, and wherein the at least one paddle marker is invisiblewithin an outline of the breast in the x-ray image; identifying at leastone region of interest in the x-ray image that corresponds to at least aportion of the at least one paddle marker; processing the at least oneregion of interest so as to derive geometric information regarding thecompressed breast image in the x-ray image; and processing the geometricinformation to derive compressed breast information.
 11. The method ofclaim 10, further comprising performing a second x-ray imagingacquisition concurrently on the breast compression paddle comprising theat least one paddle marker and the patient's breast compressed betweenthe compression paddle and the x-ray imaging receptor, wherein thesecond x-ray imaging acquisition is performed with a second x-ray dose.12. The method of claim 11, wherein the second x-ray dose is based atleast in part on the compressed breast information.
 13. The method ofclaim 11, wherein the second x-ray dose is higher than the first x-raydose.
 14. The method of claim 10, wherein the compressed breastinformation comprises breast density information.
 15. The method ofclaim 14, further comprising processing the breast density informationto derive predictor information related to cancer risk.
 16. The methodof claim 15, further comprising recording the cancer risk predictorinformation in association with image information regarding the imagedbreast.
 17. The method of claim 14, further comprising recording thebreast density information in association with image informationregarding the imaged breast.
 18. The method of claim 14, furthercomprising displaying the breast density information in association withimage information regarding the imaged breast.